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Abstract 

Two  calculated  methods  of  film  thickness  of  non 
cosine  distribution  and  two  measuring  methods  based  on 
Fabry-Perot  metal  dielectric  filters  were  presented. 
Some  calculations  and  measurements  were  shown. 

Key  Words:  Thickness  distribution,  thin  films 


0.  INTRODUCTION 


The  even  distribution  of  film  thickness  is  very 
critical  for  obtaining  thin  films  in  many  applications.  For 
example,  the  interference  filter  used  in  astronomical 
observations,  to  a  large  extent,  depends  on  the  even 
distribution  of  film  thickness  for  proper  function.  In  other 
cases,  the  unevenness  of  films  is  desirable  in  order  to 
obtain  an  appropriate  thickness  range.  Therefore, 
calculating  and  measuring  film  thickness  distribution  is  an 
essential  part  of  thin  film  technology. 

Here,  two  film  thickness  calculation  methods  will  be 
given.  In  addition,  two  measuring  methods  based  on  Fabry- 
Perot  metal  dielectric  filters  were  used  to  measure  samples 
that  have  higher  thickness  evenness  and  to  measure  samples 
that  have  great  fluctuation  in  their  film  thickness.  In  the 
end,  significant  measurements  will  be  given. 


*  Numbers  in  margins  indicate  foreign  pagination. 
Commas  in  numbers  indicate  decimals. 
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I.  CALCULATION  OF  FILM  THICKNESS  DISTRIBUTION 


Thickness  of  material  which  evaporates  from  evaporating 
source  E  and  deposits  on  a  spinning  substrate  S  can  be 
calculated  as (fig.  1): 

t=fcPV  cos^cose  ^  (1) 


In  which,  k=cm/27rcr,  a  and  m  are  the  density  and  mass 
of  the  evaporating  material  respectively,  c  is  a  constant, 
<p,  0  ,  r  and  <p  are  depicted  in  £ig.l. 


Fig.  1  The  geometric  position  of  evaporating  source  and 
substrate 


Since  cos  <p  =cos0  =  H/r, 

y*  =  H*  +Rl  +  Ll  +2RLcos<f> 


therefore  equation(l)  becomes 

*.*j; - 

ZHt+Ri  +L*+2RLcos03  * 
Assuming  p  =h/R  and  a=L/R 


p  +  3 


-d<f> 


Cl  +p2  +ol  +2ocos03 
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The  thickness  of  substrate  center  (L=0)  is 


t*  = 


kz 


r>P  +  * 


H 2  ^ 

Cl  +p2]  2 


T  =  ±  =  J1  +p 


P  +  3 


P2)  2  f* 

x  J0 


d<t> 


P  4-3 

Cl +P2  +  a2 +2iocos#J  s 


Here  T  is  film  thickness  distribution. 


If  assuming  S  =  l  +  'j°+p 


and 


Q=  2fl-. 
1  +  p2 


=— r 

JT  Jo 


d<f> 


r±3 


(2) 


CS  +  Ceos#]  * 

By  substituting  different  values  for  P,  the  film 
thickness  distribution  of  different  evaporation  sources  can 
be  calculated.  P=0  corresponds  to  a  point  evaporation 
source.  P=1  corresponds  to  Knudsen's  small  surface  source 
which  is  a  cosine  distribution.  At  the  location  where 
electron  beams  evaporate,  P  is  often  around  2-4. 

Two  simple  methods  are  used  to  calculate  thickness 
distribution  T.  The  first  is  series  expansion,  (2)  can  be 
expanded  as 


n — r - 1 


(P  +  3XP  +  5) 


S~ 


212* 


f(l)' 


(P+3)(p  +  5)(P  +  7)( P+9)  3  /  Q  \*  mm  1 

4\2*  8  \  S  /  J 


The  calculation  of  the  sum  of  the  first  8  items  is 
given  in  Table  1.  It  can  also  be  solved  as  the  following. 
Because  when  p=l,  (2)  can  be  simplified  as 


(H2  +  L2)2(H2  +  L2  +  R2) 
C(ff2  +  L2  +  R2)2-4L2R2J*/2 
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As  a  =  L/R  increases,  the  discrepancy  of  film  thickness  /401 

calculation  gets  bigger  when  choosing  the  series  expansion 

method.  Therefore,  the  second  calculation  method  —  Gauss 

integration  is  used  instead.  It  can  be  simply  described  as 

if  ,  =  (bZfl)  (b+o) 

2  2 

then  £«*<>»<“  =(nr)  2h.K«‘.» 


Here  <p  (tj_)  and  corresponding  values  of  Hj^  can  be  found  in 
relevant  mathematical  manuals. 

Suppose  0(*>  =  -^— (l  +  t)  i.e.  a=0,  b=  n 

2 

Then 

M)d<f>  =  -frj* )dt  — 2 Bif  (*«,)) 

{ml 

Therefore  equation(2)  can  be  written  as 

7"  -  i_r  _ dt _ rr  _ 1 _ 

6  LS  +  Qcos(<f>(t))2~^~  ■-«  2CS  +  Qcos(— 5-(l+t;) )]~ 

When  n=6,  the  calculated  results  are  listed  in  “Table  1. 
Apparently  it  agrees  well  with  the  results  obtained  when 
applying  Gauss  Integration. 
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TABLE  1.  THICKNESS  DISTRIBUTION  T  CALCULATED  WHEN 
p  =1  AND  P=l. 


a 

T 

m  «  s  —  1 

mmm  -2. 

0.0 

1.0000 

1.0000 

1.0000 

0.1 

l.OOSO 

1.0050 

l.OOSO 

0.2 

1.0194 

1.0194 

1.0194 

0.3  | 

1.0418 

1.0418 

1.0418 

0.4  ' 

1.0697 

1.0697 

1.0697 

0.5 

1.0991 

1.0991 

1.0992 

0.6 

1.1249 

1.1248 

1.1250 

0.7 

1.1408 

1.1405 

1.1409 

0.8 

1.1404 

1.1394 

,  1.1405 

0.9 

1.1187 

1.1164  ' 

1.1187 

1.0 

1.0733 

1.0691 

1.0731 

Key:  1-  Analysis  2-  Series  expansion  3-  Gauss  integration 


Fig. 2  depicts  the  change  of  film  thickness  with  the 
change  in  p  when  p=l.  Interestingly  enough,  the  film 
thickness  distribution  hardly  changes  when  p  is  between  1.5 
and  2.5,  therefore  it  could  be  constant  in  this  interval. 
Fig. 3  shows  the  change  of  film  thickness  with  the  change  in 


p  when  p  =1  and  3.5 


/  402 


Fig. 2  Film  thickness  distribution  at  various  p  when  p—1.0. 
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Fig.  3  Film  thickness  distribution  at  various  p  when  p-1.0 
and  3 . 5 


Sub 


Fig.  4  Metal  filters  and  their  symbols 

2.  THE  MEASUREMENT  OF  FILM  THICKNESS  DISTRIBUTION 

The  film  thickness  distribution  measurement  can  be  achieved 
by  indirectly  measuring  peak  value  wavelength  with  a  Fabry-Perot 
metal  dielectric  at  various  points. 

As  we  know,  when  (  6  +  <f>, +<pt  ) /2  =  m/T 

-^-n1d1cos0  +  01  +<f>z  =  2mx  (3) 

A 
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The  filter  has  maximum  transmittance 


Ty  =  T o/Cl  +  F-sin2  (  d  +  *»+*«_)] 
T0=T1  TJ  (1-v/^F2)S 

f=4v/r7r2/  o -✓*!**•>  > 


From  (3)  we  can  obtain 


2n1d1cos0  =  <  m-ip  )  i 


(4) 


In  which,  Tp  =  ti±^i  .  if  Ag  layer  is  symmetric,  then 


=  ♦„-! _ 2JCW, 


-&Z  =tg“ 


n2  +fc*  -  n,2 
then  (4)  can  be  simplified  as 


and  when  it  is  vertically  projected, 


2nidl  =  (m  -  tf)X0 


(5) 


If  thickness  d^  is  changed  to  d^  +  5d^  with  respect  to  the 
different  positions  of  filters,  and  the  refractive  index  between 
the  layers  ^  is  changed  to  nj_  +  6n-L  with  respect  to  different 
wavelength,  and  phase  color  scattering  is  changed  to  $  +  6i|r  , 
conseguently  the  peak  value  is  changed  to  Aq  +  6A  and  then  (5) 
can  be  written  as 

(6)  /403 

2(n,+dna  )  (dj  +  dd!  )  =  (X0+dA.) 

If  it  is  divided  by  (5),  then 


(1+*Li.)(l+r-^i-)  =  <  l  - 
n  x 


— ~~r  1  (1+-T-) 
m-w  a0 


(7) 


Let 


T  =  1  + 


from  (7), 


ddi  be  actual  film  thickness  distribution 
dj 


m-w  Ao 


T  = 


<  i+^> 


7 


(8) 


(  1  + 


<U 


)  Cl  - 


^>0 

m  -  ip 


dip  <3X  ir  1  ,  X0  dn  i  <5X 

ax  x,  JL  »,  ax  x„ 


T  =  (  1  + 


)  (  1  +  A 


ax 

X0 


in  which  a=T7^T  and  B  =  W^p~dT  correspond  to 

the  phase  change  color  scattering  and  inter-layer  refractive 
index  color  scattering. 


If  inter-layer  n-^  corresponds  to  ZnS,  A  can  be  solved 

from 


n !  =2.239 


j.  45200 
X* 


For  MgF2  underlayers,  n1=1.38,  i.e.  A=0.  B  can  be 
solved  from  the  relation  between  <p  and  i|i 


a^_  =  _j _ d± 

ax  jc  ax 


d<p  _  d<P  dn  +  30  8fc_  M  8#  dn , 

ax  8rT  ax  ax  ax  a» ,  ax 


— jK(n*  +Kin1t)?£j-+n1(ni -K*  -nt*  ) 
n^L  8X 


(n2  +  fc2  -«j2)2  +  4 K2 


in  which  the  optical  constancy  of  Ag  can  be  expressed  as 


n  =  0.525  -1 .825  xlO“8X  +  1 .75  xlO-,X2 
fc  =  3.41  +8.7  X10'*(X -546) 


(wavelength  500-700nm) 


Hence,  using  (8),  the  film  thickness  distribution  T  can 
be  easily  calculated  based  on  the  change  of  peak  value 
wavelength  with  respect  to  the  different  positions  of 
filters . 


The  method  described  above  can  also  be  used  to  deal 
with  filters  with  wide  range  of  thickness,  by  changing  the 
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wavelength  of  incoming  light  while  maintaining  the  vertical 
emittance  measurement.  If  the  evenness  of  filters  does  not 
fluctuate  a  bit,  a  more  accurate  method  is  to  use  single 
wavelength  lightening(i.e.  He-Ne  laser)  peak  value 
wavelength  which  can  be  located  by  adjusting  the  incoming 
angle. 

For  a  known  scale  m,  cos0  =  XAo >  which  is  inferred  from  /404 
(4)  and  (5).  For  a  smaller  incoming  angle, 

cos9  =  1  ~2n? 

therefore  TT'-fi 7.  <9> 

In  which  both  the  color  scattering  of  nx  and  i|r  with  respect 
to  wavelength,  and  the  change  of  Ag  layer  phase  with  respect 
to  the  change  of  incoming  angle  were  omitted.  When  taking 
these  into  consideration,  equation  (4)  becomes 

2d,(n  j  +<5n,)  (  1  -f2— *  )  *  (  m  +6K  )  ( 10  ) 


By  substituting  6n2  with  and 

oA 

dividing  by  (4),  we  find: 


0. 


d* 


_dk_ 

^0 


2n ,  *  tn-tp  dd 


de 


dtp 

'aT 


(ii) 


and 


dtp 

de 


de 


1  -A-B 

Since  5X  and  02o  have  a  linear  relationship,  therefore 
di 


m-tp  d6 
be  written  as 


de  is  a  function  of  02o  ,  i.e.  C02o,  hence  (11)  can 


0p 1  (  1  -  2«  i  ) 

K  2ni*(l  -A-B) 


(12) 


the  fraction  of  C  with  respect  to  p  and  s  is 

c  _  _c  _ _ K  (  n2  +  g2  +n  i*  ) _ 

P  ’  n  ,jr(m- tf)[4fc2n  i2 +  (n2 +K2 -n2,)2] 
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When  C=A=B=0  (12)  becomes  (9).  If  (12)  is  written  as 

<5i  -  g.2 

K  2n\  (13) 

then  n'  =n  a  [1  -  A-  B)(/l -2n  is  the  effective 

refractive  index. 

Therefore,  based  on  (13),  wavelength  change  can  be 
calculated  according  to  the  corresponding  inclination  angle 
of  the  measured  filter  maxima.  Consequently,  the  film 
thickness  distribution  can  be  calculated  from  (8). 

In  order  to  assess  the  effect  of  color  scattering  on 
the  measuring  error,  (8)  can  be  expanded  as 

T - 1  =  l-A-B  ) 

*0 


Let  f=(l-A-B).  Tm  is  the  calculated  film  thickness  that 
corresponds  to  the  measured  peak  value  drift,  i.e. 

Tm-1  =  SX/Xq 
then  T-l  =(TM-l)f 


Therefore 


=  (  JL  )  *  6tu  +  (  -U- )  ^f1 


dTt! 


3f 


ar 


=  f, 


ar  _ 

~W~Tu~1 


dr  u 

a t1  =  /*arw*  +  (Ty-  vydf1 


/405 


Assume 

Then  we  find 


Tm-1  =  0.1,  we  calculate  df2=0.014,  /*«Tjr*=7. 8xi0’# 
< Tif-1  )  *af*  =  1.4X10“ 4 
6T-  1.2X10"* 


Therefore,  the  thickness  distribution  measurement  error 
arises  from  the  second  item,  especially  the  error  in  f.  This 
shows  that  the  Fabry-Perot  metal  dielectric  filter  based 
measurement  of  film  thickness  distribution  should  take  into 
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account  the  refractive  color  scattering  among  the 
separating  layers  and  metal  layers . 

Fig. 5  are  the  set-ups  used  in  the  above  measuring 
methods . 


Fig. 5  Measuring  set-up  diagram 

Key:  1-  light  source  2-  lens  3-  sample  4-  light  filter 
5-  receptor  6-  monitor  7-  He-Ne  laser  8  self-adjustable 
parallel  light  tube  9-  sample  10-  angle  meter 
11-  receptor  12-  monitor 


3.  MEASUREMENT  RESULTS 


Using  above  methods  we  measured  3  commonly  used  film 
(ZnS,  Na3AlF6/  MgF2)  thickness  distributions.  The  following 
filters  were  used: 

AlAg8ZAglG 

AlAg8CAglG 
and  AlAg4MAglG 

In  which  Z,  C  and  M  designate  A0/4  film  of  ZnS,  Na3AlF6  and 
MgF2  • 

The  sample  is  spun  at  8  cycle/sec,  the  substrate 
temperature  is  room  temperature,  vacuum  degree  is  1-2x1 0-5 
Torr . 

Fig.  6  and  7  shows  the  measurement  at  H/R=0.9  and  1.7. 
As  we  see  from  the  figures,  the  actual  curve  diverges  from 
the  theoretical  value  to  various  extent (p=l).  At  H/R=1.7, 
the  sequence  is  Na3AlF6<ZnS<MgF2 ,  in  which  MgF2  shows  the 
greatest  diversion.  At  H/R=0.9,  the  sequence  becomes 
MgF2<  Na3AlF6<ZnS.  Different  evaporation  rate  also  affects 
film  thickness  distribution.  At  1.7,  the  sequence  is 
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ZnS<Na3AlF6<MgF2 /  in  which  MgF2  is  affected  greatly.  At  0.9, 
the  sequence  changes  to  ZnS<MgF2<Na3AlF6 .  At  high 
evaporation  rate,  the  distribution  under  these  two  given 
conditions  tends  to  approach  straight  line  T=l,  in  other 
words,  thickness  evenness  has  been  improved.  Comparing 
H/R=l . 7  and  0.9,  the  constant  value  at  H/R=1.7  is  more 
striking  than  at  0.9.  Even  comparing  the  first  time  use  and 
3  time  use  evaporation  source,  the  distribution  at  1.7  is 
more  stable,  the  sequence  is  ZnS<Na3AlF6<MgF2 ,  MgF2  is  the 
most  unstable  one,  this  is  confirmed  by  fig.  2. 

/406 

Therefore,  in  order  to  make  an  adjustable  filter,  H/R 
should  be  close  to  1.7.  In  addition,  ZnS-Na3AlF6  or  ZnS-MgF2 
film  thickness  should  be  close  .  Furthermore,  due  to  the 
great  difference  in  vapor  evaporation  between  new  and  old 
evaporation  source,  the  evaporation  source  for  low 
refractive  index  material  should  be  replaced  promptly, 
otherwise,  the  adjustable  filter  would  not  work. 
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Fig. 6  The  comparison  of  three  film  layer  thickness 
distribution  with  theory  when  H/R=1.7. 

Fig. 7  The  comparison  of  three  film  layer  thickness 
distribution  with  theory  when  H/R=0.9. 

Key:  1-  divergence. 

We  also  measured  film  thickness  distribution  of  ion 
assisted  thin  film.  We  found  that  the  evenest  film  thickness 
distribution  could  be  achieved  by  adjusting  the  distance  and 
angle  of  ion  gun  from  the  target. 

We  would  like  to  thank  Prof.  Lissberger  at  British 
Royal  University  for  his  help  in  our  study. 
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